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ABSTRACT 



A static pressure probe was tested to determine the feasibility of 
using the probe, as an integral part of a missile nose, to sense missile 
altitude. Experiments were conducted at Mach 2.0 and at Mach 1.51. 

At Mach 2.0, the static pressure probe will perform within altitude 
specifications of 25,000 feet ± 2,000 feet at angles of attack ranging 
from - 8 to + 8 degrees. At Mach 2.0, within an angle-of-attack ranging 
from 0 to 6 degrees, the probe will measure free stream static pressure 
within 4 percent; a 4 percent error in measurement is equivalent to an 
altitude error of 900 feet. The missile nose shock will remain downstream 
of the probe pressure ports for flight Mach numbers above 1.5. 
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I . INTRODUCTION 



A. BACKGROUND 

The ZEPPO Rocket is an air-launched system designed to deliver an 
expendable pulse power communications jammer at long range. To achieve 
this objective, the system is loft- launched at angles up to 45 degrees 
and at high altitude. At the proper point in the trajectory, a deploy- 
ment sequence is initiated by the fuze. At the end of this sequence, the 
expendable jammer is suspended and operated from a 17-foot diameter para- 
chute with a descent rate of 10 to 12 feet per second. The ZEPPO system 
is 5 inches in diameter and 115 inches in length, and weighs approximately 
135 pounds. The deployment system consists of: (1) a timer fuze in the 
nose; (2) an expulsion charge; (3) a drogue parachute; (4) a main para- 
chute; and (5) an expendable electronic jammer. 

The rocket motor is the production MK71 MOD 1 Zuni motor. The motor 
weighs 80 pounds, and is 73.3 inches in length and 5 inches in diameter. 

It has a thrust level of 6500 pounds force for 1.5 seconds, with a burn- 
out spin rate of 30 to 33 revolutions per second. 

The fuze, at the present time, is a modified MK375 MOD 0 timer fuze 
with a timing capability up to 140 seconds. It contains a black powder ex- 
pulsion charge. 

The drogue parachute has a cruciform geometry and is packed in a 
small bag that opens after motor separation. The main parachute is a 
lightweight 17-foot parachute adapted from the LUU-2/B flare system which 
is a standard production item. The chute is packed in a deployment sleeve 
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measuring 5 inches in diameter by eight feet in length; this sleeve dras- 
tically reduces snatch and opening loads, permitting the lightweight para- 
chute to survive high "Q" deployment without damage. 

The optimum point for deployment initiation is 25,000 feet, plus or 
minus 2,000 feet. Deployment at this altitude allows maximum payload time 
aloft, and therefore maximum operation time. System testing by the Naval 
Weapons Center, China Lake, California, has demonstrated that there is 
difficulty in attaining the optimum deployment altitude with the present 
fuzing system. 

To remedy this situation, several alternatives were considered, one 
of which was a fuzing system initiated by a pressure-sensing switch, or 
Baro-Fuze. Two alternatives were considered for measuring pressure: a 
shock-swallowing air sensor and a static-pressure probe. The Naval Post- 
graduate School was selected to evaluate the effectiveness of the static- 
pressure probe. The probe selected for this evaluation was developed by 
Pinckney [1]. 

Figure 1-1 shows the probe interior to the missile nose, where it 
remains until after launch. After rocket motor burnout, and with the mis- 
sile flying at small angles-of-attack, the probe is extended. With the 
probe extended, as shown in Figure 1-2, the aerodynamic drag is reduced; 
thus the decay of missile velocity is decreased. The range of possible 
Mach numbers was considered to be Mach 1.4 to Mach 2.0 at small angles- 
of-attack. 
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B. THESIS PURPOSE 



The purposes of this thesis are: (1) to evaluate theoretically the 
static pressure probe general designs developed by Pinckney [1] for use 
in the Mach 1.4 to Mach 2.0 range; (2) to develop pressure distributions 
along the probe, and select the probe design best suited to this Mach 
number range; (3) to evaluate experimentally the probe pressure-sensing 
performance as an integral part of the missile nose (with probe extended) ; 
and (4) to report on experimentation in the Naval Postgraduate School 
supersonic wind tunnel utilizing the existing Mach 1.4 and Mach 2.0 noz- 
zle blocks. 
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II. THEORY 



A. STATIC PRESSURE PROBE DESIGN 

Pinckney [1] discusses the design, for supersonic flow, of a short 
static pressure probe that is relatively insensitive to angle-of-attack . 

The general probe design, which is illustrated in Figure 2-1, consists of 
a conical nose followed by a parabolic tangent ogive fairing the nose into 
a second, truncated cone. The second cone mates with the cylindrical after- 
body of the probe. The parabolic formula used for the tangent ogive sec- 
tion was of the form: 



The value of Y/R^ at the junction between the tangent ogive and the 
second truncated cone was specified. The conical nose half-angle (8) and 
the second truncated cone half-angle (co) were also specified. Based on 
this information, the constants A, B, and C were calculated. The probe 
nose half-angle is dependent on the free-stream Mach number, the cone 
angle must be small enough so that shock detachment does not occur. Static 
pressure holes are located at points where static pressure equals free- 
stream static pressure. These points of pressure equality vary with the 
length and angle of the conical nose. Pinckney [1] calculated the theor- 
etical pressure distribution for this general design utilizing a blunt- 
body program coupled with a method of characteristics program. He made 
calculations for 8 = 10 degrees with u varying between 2 and 4 degrees 




n 



n 



( 2 - 1 ) 



15 



for M ot = 2.5 and 4.0 and zero angle-of-attack . Wind tunnel tests were con- 
ducted utilizing four different probe types. The first two probe types had 
3 = 10 degrees and u = 2 or 3 degrees, and were tested at M^, = 2.5 and 4. 

The second two probe types, which had 6 = 20 degrees and co = 2.0 or 3.5 
degrees, were also tested at = 2.5 and 4. Tests were conducted at attack 
angles ranging from 0 to 12 degrees. 

Figure 2-2, which is reproduced from Kuethe and Chow [2], is a graph 
of shock angle versus Mach number for various cone angles. Since the thrust 
of this research concerns Mach numbers between 1.4 and 2.0, Figure 2-2 
shows that 6 cannot be greater than approximately 15 degrees, or shock wave 
detachment will result at M^ = 1.4. Therefore, a conical nose half-angle 
of 8 = 10 degrees was chosen. Pinckney's theoretical results for 8 = 10 
degrees, which are shown in Figure 2-3, indicate that the pressure distri- 
bution which most closely approximates that of the free-stream static pres- 
sure will be obtained with oj = 3 degrees at M,* = 2.5. Pinckney's experimental 
results, which are shown in Figure 2-4, indicate close agreement with the 
theoretical calculations. Figure 2-5 illustrates Pinckney's probe design 
for 8 = 10 degrees and oj = 3 degrees. Since Pinckney's theoretical calcu- 
lations were for much higher Mach numbers, it was considered prudent, for 
research completeness, to obtain theoretical values for the static pressure 
distribution at lower Mach numbers. For consistency, attempts were made to 
obtain the computer program used by Pinckney [1] ; however, the program was 
not immediately available. Therefore, an alternative approach was sought. 

The coefficients of the parabolic tangent ogive sections for 8 = 10 degrees 
and a) = 3 degrees were readily available and were provided by Mr. Pinckney. 

The values of the coefficients are as follows: A = -0.009267312, B = 8.3676155, 
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and C = -402.88107. The coefficients were verified, and Eq. 2-1 provides 
the ogive section body shape used in the theoretical work described in the 
next section. 

B. STATIC PRESSURE DISTRIBUTION 

Gawain and Schonberger [3] discuss small perturbation supersonic flow 
theory applied to a body of revolution. The theory presented in Ref. 3 is 
essentially that of classical slender-body small -perturbation theory, with 
corrected boundary conditions in the sense that the boundary conditions 
are applied at the body surface. In more elementary theory, the boundary 
conditions are applied along the axis of the body. The use of the exact 
boundary conditions yields more accurate solutions than those obtained 
from more elementary theory. The theory of Ref. 3 was considered adequate 
for the purposes of this thesis, as only small angles-of-attack were con- 
sidered. A BASIC language computer program, suitable for use on a micro- 
processor, was developed based on the theory and calculation procedures 
in Ref. 3. The program, and a sample data run for = 2.0 at angle-of- 
attack a = 0 and 10 degrees, are given in Appendix A. 

C. RESULTS OF THEORETICAL CALCULATIONS 
1 . Theoretical Calculations 

The pressure distribution was calculated for = 1.4, 1.6, 1.8, 
2.0, and 2.5 at zero angle-of-attack (a). The M,,, = 2.5 calculations served 
as a comparison between Pinckney's theoretical and experimental results [1] 
and the theory used in this research. Figure 2-6 is a graph of P/P^ versus 
X/R n for Mo, = 2.5, g = 10 degrees, and w = 3 degrees. Comparing Figures 
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2-3 and 2-6, it can be seen that the graph shapes are the same; however, 
the theoretical results depicted in Figure 2-6 are slightly higher than 
those of Figure 2-3 (approximately 2.7 percent). Yet, in comparing Figure 
2-6 with the experimental results of Figure 2-4, only a 1 percent differ- 
ence exists . 

In addition to the calculations above, computations were made for 
angles-of-attack a = 1, 2, 4, 6, 8, and 10 degrees at M*, = 1.4, 1.6, 1.8, 
and 2.0. The calculation procedure and results will be discussed in Sec- 
tion 4, below. 

2. Probe Pressure Hole Location 

As was discussed previously, in Part A of this chapter, optimum 
probe pressure hole location depends on the free-stream Mach number. Due 
to the range of Mach numbers, a design compromise is made that minimizes 
the pressure error over the range of M OT rather than for one specific M*. 
Figure 2-7 is a graph of P/P*, versus X/R^ for M*, = 1.4, 1.6, 1.8, and 2.0, 
at zero angle-of-attack. The graph shows clearly that the minimum differ- 
ence in P/P*, for the range of M*, occurs at approximately X/R^ = 260. The 
difference in P/P*, due to variable M*, is approximately 4 percent. 

3 . Error in Altitude Due to Error in Pressure Measurement 

Pressure as a function of altitude is given by: 




where 5 = (7.6 km) 1 and P Q is the pressure at sea level. Eq. 2-2 was ob- 
tained by graphing pressure versus altitude using data from the U. S. 
Standard Atmosphere , 1976 [4], and fitting a curve to the data. Taking the 
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natural logarithm of both sides of Eq. 2-2 yields: 



In P = -E,h + tn? ^ 

Since P Q = 1 atmosphere, Eq. 2-3 becomes: 



(2-3) 



ZnP = -£h . (2-4) 

To find the error in altitude due to the difference in pressure, the deri- 
vative of Eq. 2-4 is taken, which gives: 

dP/P = -tdh. . (2-5) 

Rearranging Eq. 2-5, the altitude error becomes: 

dh = - i • dP/P . (2-6) 

With 1/C = 7.62 km and dP/P = 0.04, the value of dh is 300 meters, or 984 
feet. This error is well within the tolerance specified for probe perform- 
ance, i.e., 25,000 feet ± 2,000 feet. 

4 . Probe Pressure Readings at Angle-of-Attack 

Figure 2-8 shows a cross section of the static pressure probe. 

The static pressure hole on the windward side is number 1, while the hole 
on the leeward side is number 3. Using the principle of conservation of 
mass, the mass flow rate in, m . , must equal the mass flow rate out, m ; 
therefore: 
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»lVl * AVz * »3 V 3 A 3 ' 0 



(2-7) 



Since p 1 = p 2 = p 3 and A 1 = A 2 = A 3 , Eq. 2-7 reduces to: 

V 1 + 2V 2 + V 3 = 0 . (2-8) 

Note that Eqs . 2-7 and 2-8 require the velocity to be an algebraic quan- 
tity: a flow into the probe is defined as positive. Streeter and Wylie [5] 
demonstrate that flow through an orifice can be represented by: 

P Q - P = %CpV 2 , (2-9) 

where P Q is the pressure upstream of the orifice, P is the pressure down- 
stream of the orifice, and C is a constant for incompressible flow. Writing 
Eq. 2-9 in terms of the parameters of Figure 2-8 yields: 



P 1 


- ? 1 

P 


- ^1 V 1 C 1 




(2-10) 


P 2 


- P 1 
P ] 


- %p 2 v^c 2 




(2-11) 


P 3 


- P 1 

P 


= ^3 V I C 3 


9 


(2-12) 



where P^ is the pressure measured by the probe and C is the flow coeffici- 
ent for pressure drop across the pressure holes. Solving Eqs. 2-10, 2-11, 
and 2-12 for V 1 , V 2 , and V 3 , respectively, and substituting into Eq. 2-8, 
yields : 



[2|P 1 -P p |] 15 ± 2[2|P 2 -P ? |] % ± [2\P s -? p \] h = 0 . (2-13) 



Simplification to this form was made possible because = = P 3 anc * 

Cj = C 2 = C 3> Note that the value of C is not needed, since C cancels in 
Eq. 2-13. 
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Eq. 2-13 is nonlinear, and its eight mathematically possible cases 
are shown in Table 2.1. A plus sign (+) indicates flow into the probe; a 
minus sign (-) indicates flow out of the probe. Some cases which are mathe- 
matically possible but not physically possible can be eliminated. 



Case 

1 

2 

3 

4 

5 

6 

7 

8 



P,-P 



+ 

+ 

+ 



P -P 

2 p 



+ 

+ 




+ 

+ 

+ 



+ 



Table 2.1 Mathematically Possible Cases for Eq. 2-13 



Cases 1 and 8 can be eliminated immediately as not satisfying con- 
servation of mass. In addition, at positive angles-of-attack (i.e., probe 
tip up), the pressure P^ will be higher than P„ or P^ . Therefore, cases 2, 
6, and 7 may be ruled out as not possible. An iterative computer program 
was written to solve Eq. 2-13 for P^, given that cases 3, 4, and 5 are 
the only possible cases. A listing of the program and a sample data run 
can be found in Appendix B. 

Figure 2-9, which plots P/P OT as a function of a, shows probe theo- 
retical performance for values of a ranging between zero and ten degrees . 
The data in Figure 2-9 were calculated for = 1.4 and for M^, = 2.0. 
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III. EXPERIMENTAL APPARATUS 



A. WIND TUNNEL 

1 . Mach 2.0 Wind Tunnel Tests 

Wind tunnel tests were conducted in the Naval Postgraduate School 
4x4 inch supersonic wind tunnel using the existing Mach 2.0 nozzle blocks. 
(The Mach number of 2 is nominal; the actual Mach number is 1.92). Figures 
3-1 and 3-2 are two photographs of the general experimental apparatus. They 
present two views of the supersonic wind tunnel and its associated instru- 
mentation. Utilizing a one-fifth scale model, and regulating wind tunnel 
stagnation pressure, an actual Reynolds number of Re = 1.3 X 10 5 was at- 
tainable, which compares favorably with the Reynolds number for the full- 
scale probe at 25,000 feet (viz.. Re = 4.3 x 10 5 ) . The reference length 
for the calculation of the Reynolds number is 0.020 meters. Twenty-three 
wind tunnel tests were conducted at M*, = 1.92 with a=0, 1, 2, 4, 6, 8, 
and 10 degrees. 

Prior to any testing, the fine-mesh wire screen in the wind tun- 
nel plenum was discovered to be tom, and hanging loose in the plenum 
chamber. The screen serves to insure a flat velocity profile (flat stag- 
nation pressure profile) and to generate small-scale turbulence which dis- 
sipates prior to arrival at the test section. The plenum chamber exterior 
can be seen in Figure 3-1. Sufficient manpower to replace the plenum cham- 
ber screen was not available. Therefore, the screen was removed. To verify 
that the flow in the wind tunnel test section was uniform, an additional 
pressure tap was installed opposite one of the existing pressure taps. 
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Testing showed no significant pressure difference (approximately 0.5 per- 
cent) between the two pressure taps. Therefore, the pressure measurements 
suggest that the test section velocity profile was reasonably uniform at 
Mach 2.0. 

2. Mach 1.4 Wind Tunnel Tests 

Tests were attempted utilizing the existing Mach 1.4 nozzle blocks. 
Supersonic flow could not be established in the test section. Appendix C 
provides a sequential listing of testing attempts and corrective measures 
taken between tests, and discusses the three most probable reasons for 
testing failure. 

3 . Mach 1.51 Wind Tunnel Tests 

Since experimental data at a lower Mach number were still desired, 
various alternatives were explored. The four alternatives investigated were: 
(1) to replace the plenum chamber screen; (2) to design and manufacture 
Mach 1.5 nozzle blocks; (3) to enlarge the existing Mach 1.4 nozzle block 
second throat by machining; or (4) to attempts tests using a combination 
of Mach 1.4 and Mach 2.0 nozzle blocks. Alternatives 1 through 3 were ruled 
out due to inadequate manpower, funding, and/or research time. Alternative 
4 was attempted. Testing in this unorthodox manner was attempted in order 
to verify that the missile nose shock did not detach at a lower Mach num- 
ber. Two shock waves occur at the nose of the missile with the probe extended. 
A nearly conical shock is attached to the probe. A curved shock wave origi- 
nates due to the missile ogive; this curved, shock wave is termed the missile 
nose shock. The probe conical shock and missile nose shock merge to form 
the bow shock wave. The probe conical shock, missile nose shock, and re- 
sulting bow shock wave can be seen in Figure 3-3, which is a photograph of 
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the model mounted in the wind tunnel. The probe is oriented at zero angle- 
of-attack, and the flow is at M,,, = 1.92. Both the probe conical shock wave 
and the missile nose shock wave are attached. If the missile nose shock 
moves upstream of the probe pressure ports, erroneous pressure readings 
are obtained. 

The Mach number gradient in the wind tunnel test section resulting 
from the combination of nozzle blocks was calculated using the method of 
characteristics found in Liepmann and Roshko [6] . Figure 3-4 shows the 
Mach number gradient in the test section and the shock waves that result 
using the combination of nozzle blocks. Figure 3-5 is a photograph of the 
Mach 2.0 and Mach 1.4 nozzle blocks installed in the wind tunnel. The theo- 
retical Mach number obtained using this combination of nozzle blocks was 
M ot = 1.51 at the probe nose, M^ = 1.538 at the top, and M*, = 1.16 at the 
bottom of the wind tunnel test section. Experimental results, which are 
discussed in Chapter IV, were used to verify theory. 

B. MODEL DESIGN 

A one-fifth scale model of the static pressure probe and missile nose 
was used. Since a complete model could not be mounted in the tunnel, the 
missile nose was truncated at 40 percent of the nose length. Maximum model 
size was dictated by maximum allowed wind tunnel blockage (approximately 
6.4 percent) at Mach 1.4. The allowed wind tunnel blockage limited the 
combined frontal area of the model and mounting system. Pinckney [1] sup- 
plied the drawings for the static pressure probe. The Naval Weapons Center, 
China Lake, California, provided the formula and coefficients for the Von 
Karmen ogive missile nose. The Von Karmen ogive is given by: 
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(3-1) 



Y = -j-L • [cos _1 (l - ^0 - %sin(2cos 1 (1 - 

for the full-scale missile, where R = 2.5 inches and L = 10 inches. Figure 
3-6 is a front view of the model and double-wedge airfoil used to mount 
the model in the wind tunnel. Design compromises between the maximum mis- 
sile nose diameter and maximum wedge thickness were required to remain 
within wind tunnel blockage limits. The maximum diameter of the model nose 
used is 0.792 inches, and the maximum wedge thickness is 0.165 inches. The 
mounting wedge length is 3.900 inches, and a wedge half-angle of 9 degrees 
was utilized to prevent oblique shock detachment at Mach 1.4. The frontal 
area of the model plus mounting system was 1.005 square inches, which is 
within the allowed wind tunnel blockage. Figure 3-7 is a photograph of 
the model ogive nose and the extended static pressure probe. The ogive is 
1.200 inches in length. The pressure probe extends 0.680 inches from the 
missile nose. Figure 3-8 is a photograph of the static pressure probe 
mated to the missile nose. The four probe pressure holes are 90 degrees 
apart and 0.0135 inches in diameter, and are located 0.1325 inches from 
the probe conical nose tip. The parabolic tangent ogive and second trun- 
cated cone start at 0.0525 inch and 0.1225 inch, respectively, from the 
probe nose tip, and are made of 0.060-inch outside diameter, 0»040-inch 
inside diameter stainless steel tubing. The probe nose cone through the 
second truncated cone is made of 0.0625-inch diameter hole on centerline. 
After machining, the drill rod is mated to the stainless steel tubing 
with silver solder. Figure 3-9 is a photograph of the top view of the 
model and mounting wedge. Figure 3-10 is a photograph of the model mounted 
in the wind tunnel, and Figure 3-11 is a close-up of the model mounted in 
the wind tunnel test section. 
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C. INSTRUMENTATION 



Pressure taps in the wind tunnel wall and on the model surface were 
connected to a Giannini 12-position scanivalve with Tygon tubing. The 
scanivalve has a range from absolute zero to 100 PSIA. The scanivalve was 
connected to a Gould Statham absolute pressure transducer. Model PA732TC- 
100-350. The transducer was powered by a five-volt excitation; it has a 
pressure range from absolute zero to 100 PSIA, with a maximum frequency 
response of 8,700 Hz. A variable-gain differential amplifier conditioned 
the transducer output to an Analog Devices 4%-digit voltmeter. The voltage 
displayed is directly proportional to the pressure. Transducer calibration 
was accomplished utilizing a 60-inch mercury column and vacuum pump. For 
higher pressure, a Wallace and Tieman bourdon tube absolute pressure gauge 
was used (range absolute zero to 100 PSIA) . Schlieren photographs of the 
flow, including shock waves, were taken. A 1,000 -watt mercury vapor light 
source and a reflex-type camera converted to accept Polaroid Type 152 film 
were used. 
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IV. EXPERIMENTAL RESULTS 



A. MACH 2.0 EXPERIMENTAL RESULTS 
1. Data Reduction Technique 

Twenty-three wind tunnel tests were conducted at Mach 2.0 with 
angles-of-attack varying from - 2 to + 10 degrees. Tables 4.1 and 4.2 
depict sample wind tunnel test data for a = 0 and 2 degrees, respectively. 



Run: 4 P = 29.85 in Hg Date: 24 September 1980 
M^, = 1.92 a = 0° T = 16.2°C Time: 0902 



Reading 


Model 

(PSIG) 


Port # 1 
(PSIG) 


Plenum 

(PSIG) 


1 


-9.37 


-9.25 


+23.12 


2 


-9.37 


-9.30 


+22.82 


3 


-9.63 


-9.49 


+22.14 


4 


-9.77 


-9.77 


+19.92 


5 


-9.86 


-9.82 


+19.27 


Table 4.1 Sample Wind Tunnel Data for a = 0 


Degrees and M*, 



Run: 6 P = 


29.88 in Hg 


Date: 24 September 1980 


M ot = 1.92 


a = 2° T 


= 20.5°C 


Time: 1027 




Model 


Port n 


Plenum 


Reading 


(PSIG) 


(PSIG) 


(PSIG) 


1 


-9.12 


-9.02 


+24.75 


2 


-9.16 


-9.06 


+24.58 


3 


-9.18 


-9.10 


+24.37 


4 


-9.21 


-9.18 


+24.14 


5 


-9.35 


-9.28 . 


+23.90 


Table 4.2 Sample Wind Tunnel Data for a =2 


Degrees and M^ 
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The ratio of measured static probe pressure to wind tunnel wall 
pressure, P/P^, was calculated for each set of data points. The arithmetic 
mean, variance (N weighting) , and standard deviation (N weighting) were 
calculated using the calculation procedure found in Haber and Runyon [7] . 
Table 4.3 gives these values for various angles-of-attack, a. 



Angle-of- 

Attack 


Theoretical 

P/P» 


Experimental 

P/P- 


Experimental 

Standard 

Deviation 


0° 


1.003 


• 0.987 


±0.00861 


1° 


1.0015 


0.984 


±0.0102 


2° 


0.996 


0.987 


±0.0097 


4° 


0.976 


0.978 


±0.0122 


6° 


0.952 


0.974 


±0.008 


8° 


0.938 


0.932 


±0.013 


10° 


0.918 


0.927 


±0.016 



Table 4.3 Angles-of -Attack Versus Theoretical and 
Experimental P/P ra for = 2.0 



2. Static Pressure Probe Performance 

Figure 4-1 is a plot of the data listed in Table 4.3. Calculated 
points are indicated by the letter X, while experimental mean values are 
denoted by an asterisk (*) and standard deviations are denoted by the 
character The scale of the ordinate emphasizes greatly any deviation 

from P/Pco = 1. Note that an experimental error bar is not shown for angle- 
of-attack; such an error bar would be horizontal. Pressure measurements 
for angles-of-attack of - 1 and - 2 degrees are included in the measure- 
ments shown for a = 1 and 2 degrees . A change in a from + 2 to - 2 degrees 
had little or no effect on pressure measurements. Figure 4-3 shows that 
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probe performance is reasonably constant in the range a = 0 to 6 degrees; 
the corresponding values of mean P/P^ are 0.987 and 0.974. For angles-of- 
attack greater than 6 degrees, the measured pressure ratio is 1.6 percent 
less than that predicted by the theory. The mean value of the measured 
pressure ratio is almost exactly that of the theory at a = 4 degrees; at 
a = 6 degrees, the measured P/P^, is 2.6 percent greater than theoretical 
predictions . 

Using Eq. 2-6 with dh = 2,000 feet (609.6 meters) and 1/E, = 7.62 
kilometers, the fractional change in pressure is 0.08. Therefore, taking 
experimental uncertainty into consideration, the probe will perform within 
altitude specifications for angles-of-attack less than or equal to 8 de- 
grees at Mach 2. As stated in Chapter III, Figure 3-3 is a photograph of 
the model mounted in the wind tunnel. The probe is oriented at zero angle- 
of-attack, and the flow is at M OT = 1.92. Both the probe conical shock wave 
and the missile shock wave are visible and attached. The horizontal lines 
just above the model are due to condensation within the boundary layer on 
the inside of the test section window. The dark vertical line aft of and 
below the model is the shadow of the Tygon tubing used to measure probe 
pressure. Figure 4-2 is a photograph of the model oriented at a = + 1 
degree within the wind tunnel operating at M OT = 1.92. 

B. MACH 1.51 EXPERIMENTAL RESULTS 

Seven wind tunnel tests were conducted using the combination of Mach 
1.4 and Mach 2.0 nozzle blocks. The model and plenum pressures, as well as 
pressures at five static pressure ports located along the wind tunnel wall, 
were recorded for each wind tunnel test. The static pressure ports extended 
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from just aft of the nozzle throat to the test section. The position of 
the static pressure ports relative to the wind tunnel test section is 
shown in Figure 3-1. The static pressure measurements were used to verify 
the calculations based on the method of characteristics referred to in 
Chapter III. The measured static pressure for the seven wind tunnel tests 
was recorded, grouped by test point, and the arithmetic mean and standard 
deviation were calculated for each point. Theoretical values and actual 



experimental values were 


compared, and are 


shown in Table 


4.4. 




Theoretical 


Observed Mean 




Percent 


Test 


Pressure 


Pressure 


Theoretical 


Pressure 


Point 


(PSIA) 


(PSIA) 


Mach Number 


Difference 


P 

1 


5.97 


5.60 


1.57 


5.9 


P 

2 


6.05 


6.35 


1.56 


4.7 


P 

3 


7.41 


7.92 


1.42 


6.5 


P 

4 


8.87 


9.12 


1.29 


2.7 


P 

C 


10.00 


10.89 


1.20 


8.1 



Table 4.4 Theoretical and Observed Pressure for a 
Combination of Mach 1.4 and Mach 2.0 Nozzle Blocks 



At Mach 1.5, a pressure difference of approximately 7 percent results 
in a Mach number change of only 0.05. Therefore, the calculated Mach num- 
bers are accurate within ± 0.05, and a predicted Mach number of M„ = 1.50 
± 0.05 at the probe tip appears reasonable. Further evidence supporting a 
value of M^ = 1.5 can be seen in Figure 4-3, which is a Schlieren photo- 
graph of the model in the wind tunnel test section during one of the test 
runs. The Mach waves emanating from the top and bottom of the test section 
are oriented at 38 and 46 degrees, respectively. At Mach angle u = 46 
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degrees, the value of M*, is 1.38. The probe conical nose shock angles are 
(j> 1 = 38 (top) and <j> 2 = 46 (bottom) , indicating that the probe is at angle- 
of-attack. A rough estimate of the angle-of-attack is given by: 

*2 " *1 

a = — - ■ . (4-1) 

With 4> ^ = 38 and ij> = 46, the value of a is 4 degrees. An accurate value 
for a would require further study. An angle-of-attack between 4 and 6 de- 
grees was predicted theoretically. Since all shocks in the test section 
are curved due to Mach gradient, all Mach/shock wave angles were measured 
locally, i.e., in close proximity to the physical object generating them. 

In addition, Figure 4-3 shows a shock wave just forward of the probe tip; 
this shock wave was predicted theoretically, and is shown in Figure 3-4. 

Due to the Mach number gradient and confused flow in the test section, 
the pressure measured by the probe is not accurate, and hence is not rele- 
vant. Table 4.5, which presents sample wind tunnel test data with the 
asymmetric nozzle blocks, is included for completeness only. The unortho- 
dox wind tunnel tests did prove, however, that the missile nose shock will 
not detach at Mach numbers above 1.5. 

Figure 4-4 is a graph of distance versus theoretical Mach number and 
pressure (PSIA) in the wind tunnel test section. The theoretical pressure 
is indicated by the character "0". Pressures measured by the probe, and 
the Mach numbers calculated from the Mach waves emanating from the wind 
tunnel walls, are also indicated. Moreover, Mach waves at several points 
in the plane of the probe nose were measured and recorded on the graph. 
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Run: 4 a = 4° M* = 1.51 P = 29.98 in Hg T = 20.5°C Date: 14 October 1980 Time: 1615 
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Table 4.5 Sample Wind Tunnel Test Data for a = 4 Degrees and -1.51 



V. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

The following conclusions result from the static pressure probe study 

(1) The static pressure probe of Ref. 1, with S = 10 degrees and 
id = 3 degrees, will, as an integral part of the missile nose, 
perform within altitude specifications of 25,000 feet ± 2,000 
feet at Mach 2.0 at angles-of-attack from - 8 to + 8 degrees. 

(2) At Mach 2.0, within an angle-of-attack range from 0 to 6 degrees 
the probe will measure free-stream static pressure within 4 per- 
cent. A 4 percent error in measurement is equivalent to an alti- 
tude error of 900 feet . 

(5) The missile nose shock remains downstream of the probe pressure 
ports for flight Mach numbers above 1.5. 

(4) Since the theory predicted probe performance within 2.6 percent 
at Mach 2.0, and since the theory predicted satisfactory per- 
formance at Mach 1.4, the static pressure probe should perform 
within specifications at Mach 1.4. 

B. RECOMMENDATIONS 

The minimum operating Mach number for the static pressure probe de- 
pends on the position of the missile nose shock wave. At some lower Mach 
number, the missile nose shock wave moves upstream of the pressure ports. 
M^ is defined as the minimum Mach number for satisfactory probe operation. 



The potential use of the static pressure probe as an altitude-measuring 

device is limited by the minimum Mach number, M , and by excursions in the 

angle-of-attack . Recommended further studies are as follows: 

Cl) Replace the missing supersonic wind tunnel plenum chamber screen, 

design and manufacture a set of Mach 1.4 nozzle blocks with a 

larger second throat, and test the probe. Determine whether 

M <1.4 or M >1.4. 
m m 

(2) Once the actual flight Mach number has been determined, recalcu- 
late the theoretical distribution, if required, and reposition 
the probe pressure ports to remove the error due to Mach number 
range--an error presently inherent in probe design. 
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3° Half Angle 
Conical Section 
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Figure 1-1 Static Pressure Probe Interior to Missile Pose 
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Figure 1-2 Static Pressure Probe Both Extended and Interior to the Missile Nose 





Tangent cone section 




37 



Figure 2-1 General Static Pressure Probe Design (Reproduced from Pinckney [1] J 




Figure 2-2 Shock Angle Versus Mach Number for Various Cone Angles 
(Reproduced from Kuethe and Chow [2] J 
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Figure 2-4 Experimental Results of Static Pressure Probe at Angle-of -Attack for M ( 
0 = 10 Degrees , and <m = 3 Degrees (Reproduced from Pinckney [ 1 ] j 



2.I4D H* — 1. 33 D 




41 



Figure 2 5 Static Pressure Probe Design for $ — 10 Degrees and to = 3 Degrees (Reproduced from Pinckney [1]J 
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Figure 2-6 Theoretical Probe Static Pressure Distribution for = 2.5 , $= 10 Degrees 3 and w = 3 Degrees 
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Figure 2-7 Theoretical Probe Static Pressure Distribution for M, 

at Zero Angle -of -Attack 



Leeward Side 




Windward Side 



Figure 2-8 Static Pressure Probe Cross Section 
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Figure 2-9 Theoretical Probe Static Pressure Distribution as a Function of Angle-of -Attack (a) for M ( 
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Figure 3-1 Naval Postgraduate School Supersonic Wind Tunnel and Associated 

Instrumentation (Rear Quarter View) 
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Figure 3-3 Sehlieren Photograph of the Model Mounted in the Wind Tunnel; 

a = 0 Degrees and = 1 . 92 
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and 2.0 (Top) Nozzle Blocks 
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Figure 3-6 Front View of the Model and Double-Wedge Airfoil Used 
to Mount the Model in the Wind Tunnel 
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Figure 3-7 Photograph of the Model Ogive Nose and the Extended Static Pressure Probe 
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Figure 3-8 Photograph of the Static Pressure Pt'obe Mated to the Missile Nose 
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Figure 3-9 Top View of the Model and Mounting Wedge 
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Figure 3-10 Photograph of the Model Mounted in the Wind Tunnel 
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Figure 3-11 Close-Up View of the Model Mounted in the Wind Tunnel 
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Figure 4-1 Plot of Experimental Results for = 1.92: 
Angle-of- Attack Versus P/P w and Probe Theoretical Performance 
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Figure 4-2 Sahlieren Photograph of the Model Oriented at a - + 1 Degree In 

the Wind Tunnel Operating at = 1.92 
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Figure 4-3 Sehlieren Photograph of the Model in the Wind Tunnel Test Section: 

Mqo = 2 • 51 and a -4 Degi'ees 
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Figure 4-4 Graph of Distance Versus Theoretical Mach Number 
and Pressure (PSIA) in the Wind Tunnel 
Test Section 
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Figure A-l Sahlieren Photograph of Shook in Mind Tunnel Test Section During an Attempt 
to Obtain Supersonic Flow with the Mach 1.4 Nozzle Blocks 



APPENDIX A: PRESSURE DISTRIBUTION PROGRAM LISTING AND SAMPLE DATA RUNS 

Note: In this typed listing, and in the sample data runs, the character "°" 

signifies "greater than," and "±" signifies "less than." The ampersand (§) 
denotes exponentiation. 

00001 LPR I NT' 'SUP ERSON I C FLOW PRESSURE COEFFICIENTS ABOUT A YAWED," 

00002 t PR I NT' ' SLENDER, FOINTED, BODY OF REVOLUTION. PROGRAM BY" 

00003 LPR I NT' ' J . R . SCHONBERGER AND K.D.TILLOTSON. BASED ON " 

00004 LPR I NT' 'CALCULAT I ON PROCEDURES BY PROF. T.H.GAWAIN." 

00005 LPRINT"NAVAL POSTGRADUATE SCHOOL, MONTEREY CA. FALL 1980" 

00006 CLEAR 

00010 CLS:PI=3. 14159 

00020 DIM AC20,20),BC20,20),C(20,20),CPC20),DC20,20),FFC20),GC20),GAC2 
0),GBC20),GCC20),GNC20),GOC20),GPC20),IXC20),JOC20),JRC20),JXC20 
),RC209,RPC20),TC20,20),XC20),PRC20) 

00025 LPRINT"ENTER XI"; : INPUT XlrLPRINT XI 

00027 LPRINT"ENTER X2";:INPUT X2:LPRINT X2 

00029 LPRINT"ENTER X3"; : INPUT X3:LPRINT X3 

00031 LPRINT"ENTER Bl"; : INPUT B1:LPRINT B1 

00033 LPRINT"ENTER B2"; : INPUT B2 : LPR I NT B2 

00035 LPRINT"ENTER COEFF A";: INPUT AA: LPR I NT AA 

00037 LPRINT"ENTER COEFF B"; : INPUT BB:LPRINT BB 

00039 LPRINT"ENTER COEFF C"; : INPUT CClLPRINT CC 

00041 LPRINT"ENTER RX2"; .* INPUT RX: LPR I NT RX 

00043 LPR I NT" ENTER FREE STREAM MACH NBR"; : INPUT MrLPRINT M 

00045 MU=ATNCl/M/SQRC-l/M :: l/M+D) 

00047 RY=RX+TAN(MU) :; (X3-X2) 

00049 LPRINT"ENTER NBR OF SUB INTERVALS" INPUT N:LPRINT N:J=N-1 
00051 LPRINT"ENTER LENGTH OF SUBINTERVAL"; : INPUT H: LPR I NT H 
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00053 FOR Z=1T0J 

00055 XCZ)=Z i: H/Cl-TANCBl)/TANCMU)) 

00057 RCZ)=XCZ) :: TANCbD 
00059 RPCZ)=TANCB1) 

00061 IF X(Z)±=X1 THEN NEXT Z ELSE 63 

00063 FOR Y=ZTOJ 

00065 R=1-AA/CTAN(MU))£2 

00067 S=-1 :: C2 :: Y :! H+BB/CTANCMU))S2) 

00069 T=-1”CC/CTANCMU))S2+CY :: H)S2 
00071 X(Y)=C-S+SQRCSS2-4 :; R 5! T))/C2”R) 

00073 RCY)=SQRCAA”XCY)S2+BB”XCY)+CC) 

00075 RPCY)=C2 :; M"XCY)+BB)/C2"RCY)) 

00077 IF XCY)±=X2 THEN NEXT Y ELSE 79 
00079 FOR Z=YTOJ 

00081 X(Z)=CCRX-X2"TANCB2))/TANCMU)+Z :: H)/C1-TANCB2)/TANCMU)) 
00083 RCZ)=TANCB2)"xCZ)+RX-TAN(B2) :: X2 
00085 RPCZ)=TANCB2) 

00087 IF X(Z)±=X3 THEN NEXT Z ELSE 89 

00089 FOR Y=ZTOJ 

00091 X(Y)=RY/TAN(MU)+Y :: H 

00093 RCY)=RY 

00095 RP(Y)=0 

00097 NEXT Y 

00099 LPRINT TABC1)"J";TAB(10)"XCJ)‘;TABC20)' , RCJ)";TABC30) , 'RPCJ)" 
00101 FOR Z=1T0J 

00103 LPRINT TABC1)Z;TABC10)X(Z);TABC20)RCZ);TABC30)RPCZ) 



63 



00110 IF Z=10 THEN INPUT"PRESS ENTER TO CONTINUE";N$ 

00120 NEXT Z 
00190 B=SQRCMS2-1) 

00192 ON ERROR GOTO 1200 

00195 U=0 

00200 FOR Z=1T0J 

00205 U=U+1 

00210 FOR Y=0TOU 

00220 TCZ,Y)=SQRCCXCZ)-Y”H)S2-CB :; RCZ))S2) 

00230 NEXT Y 

00231 NEXT Z 
00235 U=0 
00240 FOR Z=1T0J 
00245 U=U+1 
00250 FOR Y=1T0U 

00260 ACZ,Y)=LOGCCCXCZ)-CY-l) :: H)+TCZ,Y-l))/CCXCZ)-Y :: H)+TCZ,Y))) 

00270 NEXT Y 

00275 NEXT Z 

00280 FOR Z=1T0J 

00290 FOR Y=1T0J 

00300 BCZ,Y)=CTCZ,Y-1)-TCZ,Y))/RCZ) 

00310 NEXT Y 

00311 NEXT Z 
00320 FOR Z=1T0J 
00330 FOR Y=1T0J 

00340 OCZ,Y)=BCZ,Y)-RPCZ) :: ACZ,Y) 
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00350 NEXT Y 

00351 NEXT Z 
00355 U=0 

00360 FPC1)=RPC1)/CC1,1) 

00370 FOR Z=2 TO J 
00380 U=U+1 
00390 FOR Y=1T0U 
00400 W=W+C(Z, Y) :: FP(Y) 

00410 NEXT Y 

00420 FPCZ)=CRPCZ)-W)/CCZ,Z) 

00430 W=0 
00440 NEXT Z 
00450 FOR Z=1T0J 
00460 FOR Y=1T0J 
00470 W=W+A(Z,Y>"FP(Y) 

00480 NEXT Y 
00490 IX(Z)=W 
00500 W =0 
00510 NEXT Z 
00520 FOR Z=1T0J 
00540 NEXT 
00550 FOR Z=1T0J 

00560 GOCZ)=(l-(l+RP(Z)S2)”(l-IX(Z))£2+(M :c IXCZ))S2) 

00565 NEXT Z 
00610 FOR Z=1 TO J 
00620 FOR Y=1 TO Z 
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00630 DCZ,Y)=CCXCZ)-CY-1)"H) :: TCZ,Y-1)-CXCZ)-Y :; H)”TCZ,Y))/C2"RCZ)&2) 

00640 NEXT Y 

00650 NEXT Z 

00660 FOR Z=1 TO J 

00670 FOR Y=1 TO Z 

00680 DCZ,Y)=-DCZ,Y)+B£2 5: ACZ,Y)+CXCZ)-CY-1)"H)"BCZ,Y)/RCZ) 

00690 NEXT Y 
00700 NEXT Z 
00710 FOR Z=1T0J 
00720 FOR Y=1T0Z 

00730 D(Z,Y)=CD(Z,Y)+B&2 :; A(Z,Y)+RP(Z>"B(Z,Y)) 

00740 NEXT Y 
00750 NEXT Z 

00760 GC0)=0:GPC1)=1/DC1,1):U=0 
00770 FOR Z=2T0J 
00780 W=0:V=0 
00790 U=U+GPCZ-1) 

00800 G(Z-1)=H”U 
00810 FOR Y=1T0Z 
00820 W=W+D(Z,Y-1) :: GP(Y-1) 

00830 V=V+BCZ,Y) :: GCY-1) 

00840 NEXT Y 

00850 GPCZ)=C1-V/RCZ)-W)/DCZ,Z) 

00860 NEXT Z 
00870 FOR Z=1T0J 
00880 U=0:V=0 
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00890 FOR Y=1T0Z 
00900 U=U+B(Z,Y) 5; GP(Y) 

00910 V=V+A(Z,Y)”GP(Y) 

00920 NEXT Y 

00930 JX(Z)=U 

00940 JRCZ)=V 5! BS2 

00950 JO(Z)=l-JR(Z)-RP(Z)”JX(Z) 

00970 NEXT Z 
00990 FOR Z=1T0J 

01010 GN(Z)=((l+RP(Z)&2) 5: (l-IX(Z))+MS2 :c IX(Z)) :: JX(Z) 

01020 GA(Z)=-1 5{ M&2”IX(Z)&2+(1+RP(Z)&2) :; (1-IX(Z))S2 

01030 GB(Z)=(M&2-1-RP(Z)S2)”JX(Z)S2 

01040 GC(Z)=(l+JO(Z))S2 

01050 NEXT Z 

01055 CLS 

01060 LPRINT"ENTER THE ANGLE OF ATTACK A (IN RADIANS)' 4 ; : INPUT ArLPRINT 
A 

01070 LPRINT"ENTER THE RADIAL ANGLE THETA (IN RADIANS)"; : INPUT T :LPRIN 
T T 

01071 LPRINT"ENTER TEMP (KELVIN)"; : INPUT TE:LPRINT TE 

01072 LPRINT"ENTER PRESSURE (N/MS2)"; : INPUT PEtLPRINT PE 

01073 LPRINT'ENTER DENSITY (KG/MS3)"; : INPUT PO:LPRINT PO 
01080 FOR Z=1T0J 

01090 CP(Z)=GO(Z)-2”GN(Z) :: SIN(A) ;: COS(A) :! COS(T)+(GA(Z)4GB(Z)”(COS(T))S2 
-GC(Z)”(SIN(T))S2)"(SIN(A))S2 

01095 PR(Z)=l+200. 941"PO :: TE :: MS2 :{ CP(Z)/PE 

01100 NEXT Z 
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01120 LPRINTTAB(l)"vJ"; TAB( 10) "X(J) M ;TABC 20) "CP(X) M ;TABC3 2) "PRESSURE RA 
TIO" 

01130 FOR Z=1T0J 

01140 LPRINTTAB(1)Z;TABC10)X(Z);TAB(20)CPCZ);TABO2)PR(Z) 

01160 NEXT Z 

01170 LPRINT"IF YOU DESIRE TO CHANGE A AND THETA ENTER Y ELSE N";:INPU 
T N$ :LPRINT N$ 

01180 IF N$="N" THEN END ELSE GOTO1055 
01200 X(Z)=X(Z)+X(Z)".001 
01210 RESUME 
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Sample Bata Run for Mach 2.0 at Angle of Attack a = 0 and 10 Degrees 



SUPERSONIC FLOW PRESSURE COEFFICIENTS ABOUT A YAWED, 
SLENDER, POINTED, BODY OF REVOLUTION. PROGRAM BY 
J.R.SCHONBERGER AND K . D . TI LLOTSON . BASED ON 
CALCULATION PROCEDURES BY PROF. T.H.GAWAIN. 

NAVAL POSTGRADUATE SCHOOL, MONTEREY CA. FALL 1980 

ENTER XI 100 

ENTER X2 260 

ENTER X3 400 

ENTER B1 .17453 

ENTER B2 .0523599 

ENTER COEFF A-9 . 26731E-03 

ENTER COEFF B 8.36762 

ENTER COEFF C-402.881 

ENTER RX2 33.856 

ENTER FREE STREAM MACH NBR 2 

ENTER NBR OF SUB INTERVALS 20 

ENTER LENGTH OF SUB INTERVAL 15 



1 


X(J) 


R(J) 


RP(J) 


1 


21.5952 


3.80776 


.176324 


2 


43.1905 


7.61551 


.176324 


3 


64.7857 


11.4233 


.176324 


4 


86.3809 


15.231 


.176324 


5 


109.828 


20.1081 


.157449 


6 


129.756 


22.9529 


.129889 


7 


148.696 


25.2279 


.111218 


8 


166.989 


27.1293 


.0971743 


9 


184.81 


28.7579 


.0859282 


10 


202.262 


30.1736 


.0765366 


11 


219.414 


31.4156 


.0684513 


12 


236.311 


32.511 


.0613283 


13 


252.989 


33.4798 


.0549372 


14 


269.503 


34.354 


.0524078 


15 


286 


35.2186 


.0524078 


16 


302.498 


36.0832 


.0524078 


17 


318.996 


36.9478 


.0524078 


18 


335.493 


37.8124 


.0524078 


19 


351.991 


38.677 


.0524078 



ENTER THE ANGLE OF ATTACK A (IN RADIANS) 0 
ENTER THE RADIAL ANGLE THETA (IN RADIANS) 0 
ENTER TEMP (KELVIN) 288.15 
ENTER PRESSURE (N/MS2) 101325 
ENTER DENSITY (KG/M&3) 1.225 
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J X(j) CP(X) PRESSURE RATIO 



1 


21.5952 


.114208 


1.31979 


2 


43.1905 


.114208 


1.31979 


3 


64.7857 


.114208 


1.31979 


4 


86.3809 


.114208 


1.31979 


5 


109.938 


.0989367 


1.27703 


6 


129.756 


.0651711 


1.18248 


7 


148.696 


.0455148 


1.12744 


8 


167.156 


.0324173 


1.09077 


9 


184.81 


.0229932 


1.06438 


10 


202.464 


.015794 


1.04422 


11 


219.414 


.0100775 


1.02822 


12 


236.311 


5.30874E-03 


1.01486 


13 


252.989 


1.29649E-03 


1.00363 


14 


269.503 


1.17926E-03 


1.0033 


15 


286 


3 . 38801E-03 


1.00949 


16 


302.498 


5.17499E-03 


1.01449 


17 


318.996 


6.62418E-03 


1.01855 


18 


335.493 


7.80804E-03 


1.02186 


19 


352.343 


8.76572E-03 


1.02454 



IF YOU DESIRE TO CHANGE A AND THETA ENTER Y ELSE NY 
ENTER THE ANGLE OF ATTACK A (IN RADIANS) .0349 
ENTER THE RADIAL ANGLE THETA (IN RADIANS) 3.1416 
ENTER TEMP (KELVIN) 288.15 
ENTER PRESSURE (N/MS2) 101325 
ENTER DENSITY (KG/MS 3) 1.225 
J X(J) CP(X) PRESSURE RATIO 



1 


21.5952 


.132234 


1.37026 


2 


43.1905 


.132234 


1.37026 


3 


64.7857 


.132234 


1.37026 


4 


86.3809 


.132234 


1.37026 


5 


109.938 


.117314 


1.32849 


6 


129.756 


.0821772 


1.2301 


7 


148.696 


.0616945 


1.17275 


8 


167.156 


.0482622 


1.13514 


9 


184.81 


.0377925 


1.10582 


10 


202.464 


.0303393 


1.08495 


11 


219.414 


.0236319 


1.06617 


12 


236.311 


.018247 


1.05109 


13 


252.989 


.0136563 


1.03824 


14 


269.503 


.0130074 


1.03642 


15 


286 


.0147623 


1.04134 


16 


302.498 


.0161278 


1.04516 


17 


318.996 


.0171903 


1.04813 


18 


335.493 


.0180226 


1.05046 


19 


352.343 


.0190157 


1.05325 



IF YOU DESIRE TO CHANGE A AND THETA ENTER Y ELSE NN 
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APPENDIX B: COMPUTER PROGRAM AND SAMPLE DATA RUN 

FOR CALCULATING PROBE PLENUM PRESSURE 



Program 

00010 CLS 

00020 INPUT "ENTER PI"; PI 

00030 INPUT "ENTER P2";P2 

00040 INPUT "ENTER P3";P3 

00050 PRINT "T=";T,"PP=";PP,"CASE=";CA 

00060 INPUT "ENTER YOUR GUESS FOR PP";PP 

00070 A=P1-PP 

00080 B=P2-PP 

00090 C=P3-PP 

00100 IF B<=0 AND C<=0 THEN 200 

00110 IF C<=0 THEN 300 

00120 IF B<=0 THEN 400 

00130 PR I NT' 'SOMETH I NG WRONG"; :STOP 

00200 T=SQR(2 5: ABS(A))-2 j: SQR(2 5C ABS(B))-SQR(2 :: ABSCO) 

00210 CA=3 

00220 GOTO 50 

00300 T=SQRC2”ABSCA))+2”SQRC2”ABSCB))-SQRC2 :: ABSCO) 
00310 CA=2 
00320 GOTO 50 

00400 T=SQR(2"ABS(A))-2 :: SQR(2”ABS(B))+SQR(2 :: ABS(C)) 

00410 CA=1 
00420 GOTO 50 



(Case 3) 
(Case 2) 
( Case 1 ) 
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Sample Data Run (Reproduced, from CRT display) 



ENTER PI ?1 .019237 
ENTER P2 ?. 96525 
ENTER P3 ?. 964412 
T= 0 PP= 0 CASE= 0 
ENTER YOUR GUESS FOR PP ?.977 
T= .174618 PP= .977 CASE= 3 
ENTER YOUR GUESS FOR PP ?.970 
T= .0131546 PP= .970 CASE= 3 

ENTER YOUR GUESS FOR PP 7.971 
T=-. 0186603 PP= .971 CASE= 3 

ENTER YOUR GUESS FOR PP ?.9705 
' T=-3.07509E-03 PP= .9705 CASE= 3 

ENTER YOUR GUESS FOR PP ?.9704 
T= 1.14933E-04 PP= .9704 CASE= 3 

ENTER YOUR GUESS FOR PP 7.970404 
T=-l . 11237E-05 PP= .970404 CASE= 3 
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APPENDIX C: SEQUENTIAL LISTING OF TESTING 



As stated in Chapter 3, Appendix C is a sequential listing of testing 
and corrective action taken between testing attempts, and a discussion of 
the three most likely reasons for the failure of the wind tunnel to start 
at Mach 1.4. Testing attempts, and the corrective actions taken between 
testing attempts, were as follows: 

(1) Test Number One was conducted with the Mach 1.4 nozzle blocks 
installed and the model mounted in the wind tunnel test section. Super- 
sonic flow could not be established in the test section. Test section 
blockage was considered to be the most probable cause. A higher percent- 
age, 28%, of test section blockage is allowed at Mach 2; therefore, the 
decision was made to complete Mach 2 testing prior to reducing model and 
mounting system frontal area. The Mach 1.4 nozzle blocks were removed. 

The wind tunnel plenum chamber screen was discovered to be tom, and was 
hanging loose in the plenum chamber. The two functions of the screen are: 
(a) to obtain a flat velocity profile (flat stagnation pressure profile) , 
and (b) to generate small-scale turbulence which dissipates prior to ar- 
rival at the test section. Sufficient manpower to replace the plenum cham- 
ber screen was not available. Therefore, the screen was removed. 

(2) The Mach 2.0 nozzle blocks were installed. To insure that the 
flow in the wind tunnel test section was uniform, an additional pressure 
tap was installed in the wind tunnel wall opposite one of the existing 
pressure taps. Testing showed no significant pressure difference (approxi- 
mately 0.S percent) between the two pressure taps. Therefore, the test 
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section velocity profile was inferred to be reasonably uniform. Twenty- 
three Mach 2.0 wind tunnel tests were conducted with supersonic flow 
achieved in the test section (25 to 30 PSIG plenum chamber pressure) . 

(3) The Mach 2.0 nozzle blocks were removed, and the Mach 1.4 noz- 
zle blocks were reinstalled. The model was mounted in the wind tunnel, 
and testing was attempted. Supersonic flow in the test section was not 
achieved. 

(4) The model frontal area was reduced by 14 percent (yielding 5.4 
percent test section blockage) by machining 0.050 inch from the top and 
bottom of the missile nose ogive. Supersonic flow still could not be 
achieved. 

(5) Model frontal area then was reduced by an additional 30 percent 
(yielding 3.5 percent test section blockage) by machining an additional 
0.060 inch from the top and bottom of the missile nose ogive and 0.015 
inch from the top and bottom of the double edge airfoil. After the reduc- 
tion in blockage, supersonic flow could not be achieved. 

(6) Four boundary- layer bleed valves were installed in both wind 
tunnel walls, just forward of the nozzle second throat, to decrease the 
required mass flow rate through the second throat. A complete description 
of the bleed valve specifications and function can be found in Habel [8] . 
Once again, supersonic flow could not be achieved. 

(7) The bleed valves and the model were removed from the wind tun- 
nel. Attempts to obtain supersonic flow in the wind tunnel test section 
failed with no wind tunnel blockage. 

(8) A fine-mesh screen was installed in the wind tunnel just for- 
ward of the nozzle first throat. With the screen, but without the model. 
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supersonic flow could not be achieved. Figure A-l is a photograph of the 
shock that was observed in the wind tunnel test section between 25 and 65 
PSIG plenum chamber pressure. 

The three most probable causes of the failure to obtain supersonic 
flow were: fa) that the model was too large, resulting in choked flow; (b) 
that the wind tunnel plenum chamber screen was missing, causing an uneven 
velocity profile; and (c) that the Mach 1.4 nozzle block second throat was 
too small, resulting in choked flow. These factors, of course, can operate 
in combination as well as singly. 

Model size was ruled out as a cause of failure due to two considera- 
tions. First, after reduction of the model's frontal area, the actual wind 
tunnel blockage, 3.5 percent, was well below the 6.4 percent theoretical 
allowed wind tunnel blockage. Second, supersonic flow could not be achieved 
even without the model in the wind tunnel. The wind tunnel plenum chamber 
screen and the nozzle block second throat, however, remained as probable 
causes of testing failure. 

The inside dimensions of the Mach 1.4 nozzle blocks were measured, 
and the areas of the first and second throats were 14.36 and 15.61 square 
inches, respectively. The ratio of second throat area to first throat area 
was thus 1.09. Using the calculation method presented in Liepman and Roshko 
[6], the theoretical minimum area ratio allowable is 1.04. Therefore, the 
Mach 1.4 nozzle block design allows a safety factor of approximately five 
percent for the boundary layer in the second throat. 

Since there is less constriction in the Mach 1.4 nozzle than in the 
Mach 2.0 nozzle, an uneven velocity profile would affect more adversely 
the flow in the Mach 1.4 nozzle. The adverse effect of the uneven velocity 
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profile may have enlarged the boundary layer in the second throat, thus 
preventing supersonic flow. 
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